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Abstract 

This note is devoted to a simple method for proving the hypocoercivity associated to a kinetic equation involving 
a linear time relaxation operator. It is based on the construction of an adapted Lyapunov functional satisfying 
a Gronwall-type inequality. The method clearly distinguishes the coercivity at microscopic level, which directly 
arises from the properties of the relaxation operator, and a spectral gap inequality at the macroscopic level for the 
spatial density, which is connected to the diffusion limit. It improves on previously known results. Our approach 
is illustrated by the linear BGK model and a relaxation operator which corresponds at macroscopic level to the 
linearized fast diffusion. 

Resume 

Hypocoercivite pour des equations cinetiques avec termes de relaxation lineaires 

Cette note est consacree a une methode simple pour demontrer I'hypocoercivite associee a une equation cinetique 
contenant un operateur de relaxation lineaire ; il s'agit de construire une fonctionnelle de Lyapunov adaptee 
verifiant une inegalite de type Gronwall. La methode distingue clairement la coercivite au niveau microscopique, 
qui provient directement des proprietes de I'operateur de relaxation, et une inegalite de trou spectral pour la 
densite spatiale, qui est reliee a la limite de diffusion. EUe ameliore les resultats anterieurs. Notre approche est 
illustree par le modele de BGK lineaire et par un operateur de relaxation qui correspond, au niveau macroscopique, 
a la diffusion rapide linearisee. 



Version frangaise abregee 

Cette note est une contribution a la theorie de V hypocoercivite, voir [10,11,8], dont le but est d'estimer 
des taux exponentiels de retour a un equilibre global pour des equations cinetiques dans lesquelles le terme 
de collision ne controle que le retour a un equilibre local. Cette question a ete partiellement resolue pour 
I'equation, hypocUiptique, de Vlasov-Fokkcr-Planck par F. Herau and F. Nier dans [6], puis par C. Villani 
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dans [11], et pour le modele de BGK lineaire, non hypoelliptique, dans [5] par F. Herau. En introduisant 

de nouveaux operateurs, nous ameliorons et nous simplifions les resultats anterieurs, tout en mettant en 
evidence la relaxation a I'echelle microscopique due au terme de collision, ici un operateur de relaxation 
lineaire, et la relaxation a I'echelle macroscopique, qui vient d'une propriete de trou spectral (voir [1]) 
d'un operateur lie a la limite de diffusion (voir [2]) et portant sur les densites spatiales uniquement. 

Soit V un potentiel exterieur sur W^, d>l. Considerons I'equation cinetique 

dtf + Tf = Lf, f = f{t,x,v),x€R^v€R^ (1) 

oil T := V ■ V X — ^ xV ■ Vv est un operateur de transport et ou I'operateur de relaxation lineaire L est 
defini par 

L/ = n/-/, Ilf:=^Fix,v),F{x,v)>0, p = p{f) := [ f dv et Pf = p{F) . 

PF J 

Sur R'' X M'' 9 {x, v), on notera dp,{x, v) = F{x, v)~^ dx dv ou F est une mesure de probabilite strictement 
positive. Norme et produit scalaire seront par defaut ceux de L'^{dp). La donnee initiale /o € L'^idp) 
est normalisee par J^j^^^^ fodxdv = 1, et la solution du probleme de Cauchy est unique car T est 
antisymetrique et L est symetrique negatif. 

Le premier cas d'application de notre methode est le cas d'un equilibre Maxwellien : F{x,v) := 
M(z;)e-^(^) avec M{v) := (27r)-'^/2 e-l"''/^, 

Theoreme 1 Supposons que V G W^^^{M.''), d > 1, est tel que : 1) j^ie~^dx = 1, 2) il existe une 
constante A > telle que /^^ \u\'^e~^dx < A J^^ \Vxu\^e~'^dx pour tout u € H^{e~'^dx) verifiant 
J^d ue~^dx = 0, 3) il existe des constantes cq > i), ci > et 9 G (0, 1) telles que AV < | |Vxl^(a;)|^ +Co 
et \ylVix)\ < ci (1 + \VxVix)\) Vx e R'i, 4) J^, \VxV\^ e'^dx < oo. 

Pour tout r] > 0, il existe une constante positive A = X{ri), explicite, pour laquelle toute solution de (1) 
dans L'^{dij) verifie : 

Vi > , \\f{t) - Ff < (1 + r?) Il/o - e-^*. 

Le deuxieme cas d'application correspond a un modele relie a Tequation de diffusion rapide avec terme 
de derive, linearisee autour de profils stationnaires de type Barenblatt. On supposera que 

/I \ -(fc+i) 

F{x,v):=uji^-\v\'' + V{x)j , V{x) = {l + \x\^f 

pour simplificr la presentation. Ici, pp = wq F''^^^'^^^, w et loq sont des constantes de normalisation, 
positives. Des choix plus gcncraux de V ainsi que les hypotheses correspondantcs seront donnes dans [3]. 

Theoreme 2 Soit d > 1, k > d/2 + 1. II existe une constante Po > 1 telle que, pour tout j3 € (min{l, (d — 
4)/(2fc — d — 2)},/?o), il existe deux constantes strictement positives C et A, explicites, pour laquelle toute 
solution de (1) dans L'^{dn) verifie : 

yt>0, \\f{t)-Ff<C\\fo-Ffe-^'. 

Sur L^{dp), a I'aidc de I'operateur de projection 11, on definit les operateurs b/ := Il{v f), a/ := b (T /), 
a / := -n (Vx/), A := (1 + a • a 0)-^ a • b et la fonctionnelle 

H{f) ■■=l\\ff + e{Af,f). 
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Si / est une solution de I'Eq. (1), alors D{f — F) := ^H{f — F) est donne par 

^(/) := (/, L /) - e (A T n /, /) - e (A T (1 - H) f,f)+e{TAf,f)+e (L /, (A + A*) /) . 

Pour simplifier les notations, on remplace / par f — F de sorte que J^^ f dxdv = 0. La preuve des 

Theoremes 1 et 2 consistc alors a montrcr que D{f) + XH{f) > 0. Elle repose principalement sur deux 
estimations. D'une part, par I'inegalite de trou spectral, le terme 

-£(ATn/,/)<-£^||n/f 

controle les termes macroscopiques a I'ordre e, c'est-a-dire proportionnels a ||n/|p. D'autre part, on 
remarque que 

2||A/f +||TA/f <||(l-n)/f. 

Lc tormc (/, L /) = — 1|(1 — 11) /|p pcrinct alors dc controlor tous les autrcs termes, et en particulier 
||(AT (1 - n))/|p que Ton evalue en considerant I'operateur dual : si (AT (1 - 11))*/ = (a • a (1 - ll))*g 
avec g = {1 + a ■ aU)^^ f, alors u := p{g)/ pf est donnee comme solution de I'equation 

P{f) = Pf « - ^ {mp V^u) (2) 

avec mp{x) := J^^ {v]"^ F(x,v) dv, et il suffit alors d'cvaluer ||(AT(1 — n))*/||^ grace a (3). Ceci rcvient 
essentiellement a etablir une estimation pour la solution de (2). Dans le cas Maxwellien, il faut 
d'abord etablir une inegalite de Poincare amelioree : il existe une constante k > telle que, pour tout 
u e H'^{e-^dx) verifiant /^^ ue-^dx = 0, k J^,i |V,iFp dx < ||V^ w||q, d'oii Ton deduit d'abord que 
Js^d |Vxl^P |Va:Wp e~^dx est borne par puis que ||V^w||o est aussi controle par Dans le cas de la 
diffusion rapide, il suffit de multiplier (2) par V^~^/^u et par V Au, puis d'effectuer quelques integrations 
par parties, pour controler en definitive ||(AT (1 — n))*/||^ par ||/||^, ce qui permet de conclure. 



1. Introduction 



A fundamental question which goes back to the early days of kinetic theory is to estimate the rate of 
relaxation of the solutions towards a global equilibrium. This is not an easy issue since the collision term 
responsible for the relaxation acts, in most of the cases, only on the velocity space. Rates of convergence 
have been investigated in many papers for the so-called homogeneous kinetic equations, but understand- 
ing how the transport operator interacts with collisions to produce a global relaxation is a different and 
much more recent story. The point is to understand how the spatial density evolves towards a density 
corresponding to a distribution function which is simultaneously in the kernels of the collision and trans- 
port operators, a property of the stationary solutions of many kinetic equations. There is an obvious link 
with diffusion or hydrodynamic limits. A key feature of our approach is that it clearly distinguishes the 
mechanisms of relaxation at microscopic level (convergence towards a local equilibrium, in velocity space) 
and macroscopic level (convergence of the spatial density to a steady state), where the rate is given by 
a spectral gap which has to do with the underlying diffusion equation for the spatial density. See [3] for 
more details. 

First non constructive results were obtained by Ukai et al., see for instance [9]. Constructive methods 
inspired from hypoelliptic theory (see e.g. [7]) were then brought into the field of kinetic theory by F. 
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Herau and F. Nier, see for instance [6] in case of the Vlasov-Fokker-Planck equation. In a recent paper, 
[5], F. Herau studied with such tools the case of an operator of zcroth order in the derivatives, which is 
known in the kinetic Uterature as the hnear Boltzmann relaxation operator. Our approach is done in the 
spirit of [5] but in a simplified framework for which the order of the operator plays no role. Moreover, 
explicit estimates on the relaxation rate easily follow, weaker assumption on the external potential than 
in [5] are needed and the method applies to more general relaxation operators of which we shall give an 
example. This example is based on kinetic equations which have been studied in [2] and give equations 
of fast diffusion in the diffusion limit. 

The hypoelliptic theory is mainly focused on the regularization properties of the evolution operator, 
but in some cases the hypoelliptic estimates also imply a result of relaxation to equilibrium. However 
both questions are independent and have to be distinguished. In the hypocoercivity approach, the purpose 
is centered on the asymptotic behavior and the quantification of the relaxation rates. More precisely our 
goal is to construct a Lyapunov functional, or generalized entropy, and establish an inequality relating 
the entropy and its time derivative along the flow defined by the evolution equation. To establish the 
inequality is then equivalent to prove an exponential rate. Such an approach has systematically been 
tackled by C. Villani, see [10,11], and has been successfully applied to various models, see for instance 
[8]. It is also related to recent works on non-linear Boltzmann and Landau equations, see 6.17. [4]. In this 
note, we develop a new approach based on operators with less algebraic properties than the ones of the 
hypoelliptic theory, but which arc better adapted to the micro-macro decomposition of the distribution 
function and give a much simpler insight of the mechanisms responsible of the relaxation at both levels. 
We illustrate our approach on two examples: the linear BGK and the linearized fast diffusion models. We 
refer the interested reader to a forthcoming paper, [3], in which the theory will be developed at a more 
general and abstract level. 



2. Main results 

Let y be a given external potential on M*^, d> 1, and consider the kinetic equation 

dtf + Jf=lf, / = /(i,x,z;), xGR^ wGK^ (1) 
where T := v ■ Vx — '^xV • Vt, is a transport operator, and the linear relaxation operator L is defined by 

L/ = n/-/, Ilf~^F{x,v), p = p{f):= [ fdv and Pf = p{F) 
Pf J 

for some function F{x,v) > which only depends on \v\'^ /2 + V{x). On R'' x R'' 9 ix,v), we consider 
the measure dp{x,v) = F{x,v)~^ dxdv where F is a positive probability measure. Unless it is explicitly 
specified, the scalar product and the norm are the ones of L^{dp): {f,g) = /udxR-^ f S ^t^- and ||/|p = 
(/, /). Throughout this paper, Eq. (1) is supplemented with a nonnegativc initial datmn /o G LF'{dp) such 
that /jjdxRti fa dxdv = 1. We shall assume that a unique solution globally exists. This is granted under 
additional technical assumptions, see for instance [2]. The goal of this note is to state hypocoercivity 
results in the two following cases. 

2.1. Maxwellian case 

We assume that F{x, v) := M{v) e'^^^^ with M{v) := (2tt)-''/^ e-l''l'/2, where V{x) = Cfc (1 + \x\^)''/'^ 
for > 1 and Ck is appropriately chosen or, more generally, satisfies the following assumptions: 
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(HI) Regularity: V G Wi^^iR^). 
(H2) Normalization: J^^ e~^dx = 1. 

(H3) Spectral gap condition: there exists a positive constant A such that |up e~^dx < A J^^ | V^wp e"^dx 

for any u € H^{e~^dx) such that ue~^dx = 0. 
(H4) Pointwisc condition 1: there exists co > and 6 G (0, 1) such that AV < | |Va;V(x)p + cq Va; e M*^. 
(H5) Pointwise condition 2: there exists ci > such that \\/lV{x)\ < ci (1 + |Vx^(a;)|) G R''. 
(H6) Growti condition: J^^^^ | Va;Fp e~^da; < oo. 

Theorem 2.1 For ony ?? > 0, there exists an explicit, positive constant A = A(?7) such that, under the 

above assumptions, the solution of (1) satisfies: 

Vi>0, \\f{t)-Ff<{l + v)\\fo-Fre->^*. 

2.2. Fast diffusion case 

For some /3 > to be specified later, we assume that 

F{x,v):=u;i^-\v\^ + V{x)j , V (x) = {l + \x\^f 

where a; is a normahzation constant chosen such that /gd^Rd Fdxdv = 1 and pp = ujq V'^/^~'^~^ for some 
Wo > 0. More general choices for V and corresponding assumptions can be found in [3]. 

Theorem 2.2 Let d>l,k> d/2 + 1. There exists a constant /3o > 1 such that, for any /? G (min{l, {d — 
4)/(2fc — d — 2)},/3o), there are two positive, explicit constants C and A for which the solution of (1) 
satisfies: 

> , \\f{t) - Ff < C ll/o - Ff e-^*. 

2.3. A Lyapunov functional 

On L'^(dfi), n is the orthogonal projection onto the space of local equilibria. Let us define 
bf:=U{vf), a/:=b(T/), a/:=-n(V,/) and A := (1 + a • a H)"! a • b . 
In the definition of A, we take the product coordinate by coordinate. These operators can be rewritten as 

bf=— fvfdv, af = —lv,-[v(i^vfdv + p{f)W,v], and a/ = -— V,p(/). 
Pf J pF \ J / Pf 

Let us also note that AT = (l + a-aH)"! a-a and a H/ = j^^^x(^^^ where := J^^ \v\^F{-,v)dv 
= d Jjjd \vi\^ F{-, v) dv for any i = 1, 2. . .d. Define the functional 

H{f) :=i||/f + £(A/,/). 
The operator T is skew-symmetric on L'^{dii). If / is a solution of Eq. (1), then 

±H{f-F)=D{f-F) 

with := (/, L /) - £ (A T n /, /) - £ (A T (1 - n) /, /) + e (T A /, /) + £ (L /, (A + A* ) /) . 
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The proof of Theorems 2.1 and 2.2 entirely relies on the following estimate with r] = 2e/{l — e). 

Proposition 3 Under the assumptions of Theorem 2.1 or 2.2, for any e > small enough, there exists 
an explicit constant A = A(£) > such that D{f — F) + \H{f — F) <Q and liminfe^o > 0- 



3. Proofs of Proposition 3 

To simplify the computations, we replace / by / — F. Therefore, from now on we assume that = 
/M<ixR<i fdxdv = {f,F). By definition of 11, /^^(n/ - f)dv = 0. We have obviously {Lf,f) < -||(1 - 
n) / IP, and using the identity exy < | a;^ + |^ t/^, we get the estimates 
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-e (AT (1 - n) /, /) = -£ (AT (1 - n) /, n/) < I HAT (1 - n) /f + ^||n/f , 

^ Z C2 

£ (T A /, /) = £ (T A /, (1 - n) /) < I ||T A /f + 1 11(1 - n) /f , 
£ ((A + A*) L /, /) < £ 11(1 - n) /f + £ II A /||^ 

for any C2 > 0. Here we have used the identities AT (1 - 11) = 11 AT (1 - IT) and T A = (1 - 11) T A, 
which are respectively consequences of the fact that the range of A is contained in 11 {dji) and that 
TA/ = vFVx{p{^f)/pF)- We moreover observe that (a • aH/,/) = ^ J^^ \Vx{^)fmFdx. Let g = 
A/, u = p{g)/pF- An elementary computation shows that —VxJ^dVfdv = ppu — ^ V^; (m^ Vc^u), 
from which one can deduce, after a few steps that we shall omit here, that HA/p = J^^ \u\'^ ppdx and 
||TA/|p = 2 /gd IVxul^mpdx are such that 

2||A/f +||TA/f <||(l-n)/f. 

We have proved that 

^(/)<-(i-^£) ||(i-n)/|p-£(ATn/,/) + |||AT(i-n)/f + ii||n/|p. 

To complete the proof of Proposition 3, it remains to estimate from above — (ATII/, /) and ||AT(1 — 
n)/||^. As for the second of these two terms, we actually estimate ||(AT(1 — n))*/|p as follows. Using 
(AT(1 -n))*/ = (a • a(l -n))*c/ with g = (1 + a •an)"V, we first observe that 

P{f) = Pf w - ^ Vx {mp Vxu) (2) 

where u = p{g)/pF- Let qp := J^^ \vi\'^ Fdv, Uij := d'^u/dxidxj. After some elementary but tedious 
computations, we also get 



||(AT(i-n))*/f = ^ E I 
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(2 dij + l)qF - ^ij ] "^ii +2(1- qf m-,- 



dx (3) 



Case of Theorem 2.1. In the Maxwellian case, various simplifications occur. With pF = e ^ = 2 ''^f ~ 1F' 
(2) becomes 

p(/)=we-^-V,(e-^V,u) (4) 
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and it follows from (H3) that 

(ATn/,/)>^||n/f. 

On the other hand, from the above computation, 

||(AT(l-n))*/f <2 ^ I \u,jf e-'^dx. 

Let IImIIo := J^d \u\^e~^dx and W := IV^l^l. By multiplying (4) by u, we get after an integration by 

parts that \\u\\q + ||Vj;iiI|Q < l|n/p. By expanding the square in | V2;('U e~^/^)p, one can prove using 
(H3) and (H4) that the following improved Poincare inequality holds, with k = (1 — 0)/(2 (2 + Acq)): 

K\\Wu\\l<\\V.u\\l (5) 

for any u € H^{e~^ dx) such that J^^^ ue~^dx — 0. 

Multiply (4) by u and integrate by parts. By (H5), we get 

\\Wu\\l + \\WWM\l - 2ci (||V.t.||o + \\WVM\^) ■ \\Wu\\o < J + I W^ff ■ (6) 

Applying (5) to Wu — jg<j W ue~^dx, we get 

k||W^^u||o< j \yx{Wu)\^ e-^dx + 2K j Wue'^dx j W^ue'^dx. 

On the one hand, by the Cauchy-Schwarz inequality, f^^ W ue~^dx < \\W\\q \\u\\q =: a, and on the other 
hand, W^ue'^dx < a ||M^||o+ IIW'^ullo' so that 2 J^^Wue-^dx 4^ W^ue'^dx can be bounded 
by i ||VK2 y||2 ^_ 2 ||vt^||4 ||u||2. Notice that \\W\\o is bounded by (H6). As for the other term of the r.h.s., 
we can simply write that J^^ |Va;(W^'u)p e^^dx is bounded by 2 Va;M||o +4cf (||u||o + ||W^u||o) using 
(H5). Hence we have 

K WW u\\l < 4 \\Wy,u\\l + 8cf (Ml + \\Wu\\l) +4K\\Wro Ml ■ 

Combined with (6), this proves that, for some C3 > 0, 

||w^v,u||o<c3||n/||. 

By multiplying (4) by Au and integrating by parts, we get 

||V2^.||2-(pFV,«||o+||^/||)||V2^i||o< ||T4-u||o||V,^.||o. 

Altogether, this proves that 1 1 ( A T ( 1 — H) ) * / 1 ^ < C4 1 1 / 1 ^ for some C4 > and, as a consequence, 1 1 ( A T ( 1 — 
n))/|p < C4 ll/f . Since (1 - n)^ = 1 - IL, we finally obtain 

||(AT(i-n))/f <C4|| (i-n)/f. 

Summarizing, with Ai = 1 — ^ C2 C4 — 5 e/2 and A2 = — we have proved that 

i)(/)<-Ai||(i-n)/|p-A2||n/f. 

With C2 = ae, a > i(l + l/A) and £ > small enough, Ai and A2 are positive and the result holds with 
A = min{Ai, A2}. The explicit expression of C4 can easily be retraced in the above computations. 



Case of Theorem 2.2. In the fast diffusion case, we only sketch the main steps of the proof. For p = 0, 1, 2, 

let wl := LOa T/p-?, where q = fc + l-d/2, := pp and V{x) = (l + \x\'^Y . Define \\u\\1 = /^^^ |u|2w2rfa;. 
Notice that pp € ^^(IR'') means I3{d - 2k -2) + d < 0. This is the case if /?(c; + 2 - 2A:) + d - 4 < and 
/? > 1. The proof in the Maxwellian case can be adapted as follows. Eq. (2) can be rewritten as 

2 

p{f) =wlu- {wl V^u) (7) 

and (H3) is replaced by the Hardy-Poincarc inequality, sec [1]. l|?il|o < A HVaiiHf for some A > 0, imder 
the condition J^^uw^dx = 0. This holds true if > 1. The fact that (ATH/,/) > ||n/|p then 
follows. We also need the following Hardy-Poincare inequality 



/ 



f L,v"+^-''-^ udxY 1 

ya+l-q-i I |2 , _ J_ < i / ya+l-q |y |2 , 



which is responsible for the condition P < Po{5), S > 0. Observe that (7) multiplied by u gives, after an 
integration by parts, ||u||q + {q — ||Va,M||i < ||n/|p. By multiplying (7) by u with a := 1 — 1/(3 
or by y Au and integrating by parts, we find directly that ||V^u||2 is bounded by ||n/||^. Computations 
which are quite similar to the ones of the Maxwellian case then allow to conclude. More details will be 
given in [3]. 
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